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ABSTRACT 



% 



Numerical models of disc galaxy formation predict the existence of extended, hot 
(T ~ lO^K) gas haloes around present day spirals. The X-ray luminosity of these 
haloes is predicted to increase strongly with galaxy mass. However, searches for their 
X-ray emission have not been successful so far. 

We calculate the all sky O VII column density distributions for the haloes of 
three Milky Way like disc galaxies, resulting from cosmological high-resolution, N- 
body/gasdynamical simulations. We perform calculations both including the disc gas 
and without it, so the disc contribution to the column density is quantified. It is 
found that the column densities estimated for Milky Way-like galaxies arc just below 
the observational upper limit, making a test of the hot halo paradigm likely within 
observational reach. 

Key words: galaxies: formation - galaxies:haloes - galaxies: spiral - X-rays: galaxies. 



1 INTRODUCTION 

Numerical simulations of disc galaxy formation predict that 
disc galaxies, not only in groups or clusters, but also in isola- 
tion, should reside in extended, hot gaseous haloes {T ~ 10® 
K). This gas is accreted from the intergalactic medium and 
cools out as it settles onto the disc, thus contributing to the 
formation of the galaxy. The mass of the predicted halo gas 
is comparable to the total baryonic mass of the galactic disc 
and coul d be the answer to th e 'missing galactic baryons' 
problem (|Sommer-Larsenll2006l '). 

Gas at temperatures T ~ 10® K is best studied in the 
X-ray regime. The X-ray luminosities of the galactic haloes 
are predicted to depend strongly on the characteristic cir- 
cular velocity of the galaxy, Lx ex: V^ (JToft et al.ll2002l ) and 
should be detectable for the most massive galaxies. How- 
ever, searches for this X-ray emitting gas have so far been 
unsuccessful for quiescent disc galaxies. Given the Lx oc Vj 
sca ling, the most restrictiv e test has been been performed 
bv iRasmussen et al.l ((20091), who observed the halo of the 
Vc ~ 307 km/s, quiescent, edge-on disc galaxy NGC 5746, 
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and obtained an upper limit compatible with, but not far 
from, theoretical predictions. 

Plenty of evidence of a hot dilute halo around the 
Milky Way exists at present: 1) the head-tail structure 
of many of the high-velocity HI clouds (HVCs) indicates 
that they are moving i n a medium much more hot and 
dilut e than themselves (jBriins et al.l l2000l : iGangulv et al.l 
120051), 2) the absence of HI in Local Group dwarf galax- 
ies within ~270 kpc distance of e ither the Milky Way 



or M31 jGrcevich fc Putrnanl |2009| ). 3) the Magellanic 



Stream (IStanimirovic et al 



of the LMC (|Mastropietrd 12003 ') an d 4) high-velocity O 



2008 ) , given the revised orbit 



VI line absorption in quasar spectra (jSembach et al.l 120031 : 
IWakker et al.ll2005l ). O VI probes gas at temperatures near 
10^ K, so high-velocity O VI associated with hydrogen in- 
dicates a thermally unstable interface between the high ve- 
locity cool clouds and a hot surrounding medium. More- 
over, cool-out of hot halo gas, at a rate of ~ 1 Mq/jt is 
indicated by Silll absorption by high- and intermediate - 
velocity clouds in the spectra of AGNs (|Shull et aLlbOOgJ ). 
Finally, the warm-h ot phase of the intergalactic medium 
(jNicastro et al.ll2005l ) provides a reservoir of hot and dilute 
gas at galactic distances. 

Six times ionized oxygen is indicative of temperatures 
around 10® K so, considering the high cosmic oxygen abun- 
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Table 1. Numerical and physical characteristics of the galaxy simulations at 2=0: mass of DM/(high-res)gas/star particles and the 
respective gravitational (spline) softening lengths; total number of particles, total number of SPH particles and total number of SPH 
particles after removal of disc etc.; virial mass and radius of the galaxy halo; characteristic circular speed of the galaxy 
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dance, the O VII line is the best way to check for the exis- 
tence of a hot gas halo around the Milky Way on the basis of 
absorption line studies. Absorption in X-ray s has been rou- 
tinely observed at zero redsh i ft for some time (INicastro et al.l 
l200alRasmussen et al.ll2003l : lFang et al.ll2006l '). but the ques- 
tion of where exactly this gas is lo cated remains under study. 
iBregman fc Llovd-Davied (120071 ) compared X-ray absorp- 
tion between lines of sight in different directions on the sky, 
and concl uded a galactic h alo origin for the absorbing gas is 
plausible. I Yao et al.l (|2008l ) made a search for the X-ray halo 
of the Milky Way by comparing O VII line absorption be- 
tween spectra of objects at different distances and positions 
relative to the galactic plane. The difference between absorp- 
tion in sightlines that mostly probe the disc and others that 
also probe the halo provides an upper limit on the contribu- 
tion of the halo gas to the total absorption. From this study, 
an upper limit of N{OVII) < 5 • 10^^ cm"^ was found for 
the sig htline considered through the galactic halo. lYao et al.l 
l|201[| ) calculated a somewhat lower column-density contri- 
bution of the putative galactic haloes to X-ray abroption by 
stacking observations towards AGNs and using the interven- 
ing galaxies as tracers. Their study gave an upper limit to 
the ionic column densities coming from galaxy haloes rang- 
ing from Zog(Af (OF//)) < 14.2-14.8, but we note that these 
constraints are less restrictive, since the sight-lines will, in 
general, pass through the haloes at considerable distances 
from the galaxy centres - see also the discussion in section 
4. 



In this paper we show that there is no discrepancy be- 
tween these observational limits and predictions from galaxy 
formation simulations. Given information on gas density, 
temperature, oxygen abundance and UVB radiation field 
from the s imulations, we use the (photo-) ionization code 
CLOUDY (JFerland et al.ll 19981 ) to determine the O VII col- 
umn density along all lines-of-sight (from positions in the 
disc corresponding to the location of the sun in the Galactic 
disc) in the simulated galaxies. We perform the analysis with 
and without the contribution from the disc gas included. 



The structure of this paper is as follows: In section [2] we 
describe the method we used for our calculations, in section 
[3] we present the results, and in section |4] we discuss the 
results obtained. 



2 NUMERICAL METHOD 

2.1 The N-body/gasdynartiical galaxy formation 
simulations 

For our analysis we used three simulated. Milky Way like 
disc galaxies, resulting from cosmological galaxy formation 
and evolution simulations based on a flat A.CDM cosmology 
{^m =0.3,nA =0.7). 

The code used for the simulations was a signifi- 
cantly improved version of the TreeSPH code, which has 
been used previously for galaxy formation simulations 
(|Sommer-Larsen et al.l (|2003t ). in the following SLOP). The 
main improvements over the previous version are: (1) The 
"c onservative" entropy equatio n solving scheme suggested 
bv lSpringel fc Hernguistl (|2002l ) has been adopted. (2) Non- 
instantaneous gas recycling and chemical evolution, tracing 
10 elements (H, He, C, N, O, Mg, Si, S, Ca and Fe), has 
been incorporated in the code following Lia et al. (2002a,b); 
the algorithm includes supernovas of type II and type la, 
and mass loss from stars of all masses. (3) Atomic radia- 
tive cooling depending both on the metal abundance of 
the gas and on the meta- galactic UV field, modelled af- 
ter |Haardt_&.^Iadau| ([199y) is invoked, as well as simplified 
treatment of radiative transfer, switching off the UV field 
where the gas becomes optically thick to Lyman limit pho- 
tons on scales of -^ 1 kpc. 

Simulations of three disc galaxies, of mass and size com- 
parable to that of the Milky Way, were used for the present 
work. The galaxies selected, denoted K33, K26 and K15 in 
the following, represent "field" galaxies (SLGP), and have, 
at 2=0, characteristic circular speeds of Vc = 180, 207 and 
245 km/s, and virial masses of 0.4-0.9x1012Mq. 

The galaxies (galaxy DM haloes) were initially selected 
from a cosmological, DM-only simulation of box-length 10 
h^^ Mpc (comoving), and starting redshift Zi=39. Mass and 
force resolution was increased in Lagrangian regions enclos- 
ing the galaxies, and in these regions all DM particles were 
split into a DM particle and a gas (SPH) particle accord- 
ing to an adopted universal baryon fraction of /b=0.15, in 
line with recent estimates. Particle numbers for the three 
TreeSPH simulations were in the range 0.3-1.1x10^ (com- 
parison to a lower resolution simulation consisting of 1.5x10^ 
particles will be briefly made). 

For the galaxy of Vc = 180 km/s, mgas=7Ti,= 9.3x10* 
and mDM=5.2xlO^ /i^^Mq. Moreover, gravitational (spline) 
softening lengths of egas=e. = 190 and eDM=340 /i^^pc, re- 
spectively, were adopted. For the galaxies of Vc = 207 and 
245 km/s, mgas=m*= 7.3x10^ and mDM=4.2xlO'^ /i'^Mq, 
and egas=e.=380 and eDM=680 h~^x>c. The gravity (spline) 
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softening lengths were fixed in physical coordinates from 
z=6 to 2=0, and in comoving coordinates at earlier times. 

A Kroupa IMF was used in the simulations, and 
early rapid and self-propagating star-formation (sometimes 
dubbed "positive feedback") was invoked (SLGP). Finally, 
in order to enable some reuse of previous work, values of 
h=0.65 and as=1.0 were employed in the cosmological sim- 
ulations. Test simulations adopting h = 0.7 and as — 0.9 
give very similar results, sin ce they produce in to tal only 
about 5 — 10% more hot gas (|Sommer-Larsenll2006h . 

In order to increase the resolution of the (dilute) halo 
gas, and enable the formation of a halo multi-phase medium 
(though still dominated by the hot gas phase; Sommer- 
Larsen 2006), for galaxies K15 and K26 (which were run 
at eight times lower mass resolution and two times lower 
force resolution than K33), each gas particle of T > 20000 
K, was split into 8 particles 0.7 Gyr before the end of the 
simulation (i.e., z=0). The simulations were then continued 
for 0.2 Gyr, and, for galaxies K26 and K33, each gas parti- 
cle of r > 20000 K was then again split into 8 gas particles. 
The resulting ultra-high resolution gas particle mass was 
1.2x10* h~^MQ, and gravity softening length egas=95/i~^pc 
for the two simulations, which were then continued for an 
additional 0.5 Gyr, till z=0. This allows the development 
of a multi-phase halo gas structure, including neutral gas 
clouds ( "high- velocity clou ds") — mo r e det ail is given in 
ISommer-LarsenI (|2006l ) and lPeek et"aLl (|2008i '). 

Finally, to test for effects of numerical resolution, a sim- 
ulation of galaxy K33, but at the (lower) resolution of galax- 
ies K26 and K15, was also run. Moreover, as detailed above, 
gas particles of T > 20000 K were split into 8 particles 0.7 
Gyr before the end of the simulation and yet again split into 
8 particles 0.5 Gyr before z—0. 



Figure 1. O VII fractions, /y, calculated using CLOYDY and 
invoking a z=0 Haardt & Madau background UV field. The efli'ect 
of the UV field becomes significant at hydrogen densities nu J; 
10"'' cm-3 
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2.3 Calculating O VII column densities 

For our analysis we used SPH simulation data for three spiral 
galaxies. We calculated column densities as, 



N{OVII) = J2'^ovii-d , 
with 



(1) 



2.2 Oxygen fractions from CLOUDY 

We calculated the fraction of oxygen in each i onisation state 
using CLOUDY version 7.2, last described bv lFerland et al.l 
([1993). CLOUDY is a spectral synthesis code, designed to 
simulate astronomical plasmas. 

We used the code for a plane-parallel geometry setup, 
assuming constant temperature o f the gas and including a 
background UV ionizing field from lHaardt fc Madaul ([1997]), 
adapted to z = as described in CLOUDY's user guide. The 
calculation was done for a single zone and iteration was al- 
lowed to continue until convergence. We thus created a look- 
up table of O VII fractions for hydrogen densities ranging 
from logn_ff — —10 to log7i_ff = —3 and gas temperatures 
from logT = 3.5 to logT = 7. The gas was assumed to be 
optically thin to the UV radiation. Self shielding of the cold 
and dense gas was approximated by "turning ofF the UV 
background for temperatures T < 20000 K and densities 
uh > W~^ cm~^ — this condition, at z ~ 0, effectively 
mimics the (simplified) UVB radiative transfer invoked in 
the galaxy formation simulations. 

Without a background radiation field and for thermo- 
dynamical equilibrium, ionisation fractions only depe nd on 
the g as temperature, according to the Saha equation (jSahal 
Il92ll 'l. However, with the inclusion of the UV field, the ion- 
ization state also depends on the density. The dependence 
of the O VII fraction on density and temperature, at z—Q, 
is shown in Figure [l] 
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Zo,G ■ fr 
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where uh is the hydrogen number density in cm~^, Zo,sim 
is the oxygen abundance in uni ts of solar, Zo,p , is th e solar 
oxygen number abundance from lAsplund et al.l ((20051), d the 
distance in cm between two points on a line of sight, /r is the 
fraction of oxygen in the sixth ionisation state and the sum is 
over a set of equidistant points along a given line of sight, nn 
and Zo,si?-n were given by the simulations, while the oxygen 
fractions were calculated using CLOUDY, as described in 
section 12.21 

For each of the simulated galaxies we drew lines of sight 
from within the disc into the halo, covering all the sky. The 
starting point of each line was placed on the midplane of the 
disc, approximately 8 kpc from the centre and the integra- 
tion was carried into the halo, up to a distance of 500kpc. 
The choice of 500kpc is discussed in section [S] 

Each line of sight was represented by 150 equidistant 
points and to each point p, a smoothing length was assigned, 



/i^= -(dso + d5i) 



(3) 



where hp is the smoothing length of the line point and dso 
and dsi are the distances of the 50th and 51st closest SPH 
particles to the point, respectively. This choice was made 
for consistence with the SPH calculations, which use the 50 
closest neighbours of a particle to evaluate interaction forces. 
All quantities of interest were then averaged over a sphere 
of radius 2/ip-, centred on the point p, using a cubic spline 
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weighting function that mimics the SPH kernel function, and 
attributed to the point. 

Specifically, for each SPH particle, on the basis of the 
density and temperature, the OVlf fraction, f-j, was calcu- 
lated using the CLOUDY based lookup table we have cre- 
ated. Using the oxygen abundance, novii was then deter- 
mined. These values were then interpolated onto the indi- 
vidual line-of-sight points, as described above. Finally, eq[T] 
was used to calculate N(OVII) along a given line-of-sight. 

For choosing the number of points to represent a line of 
sight, we compared the results of 20000 line of sight calcu- 
lations, each done using 150 points and then repeated using 
500 points and found only insignificant differences in the 
column densities. Higher line-of-sight resolution than this is 
not justified, as we would then interpolate particle proper- 
ties over distances smaller than the SPH smoothing lengths. 



3 RESULTS 

Not only hot halo gas may contribute to the O VII col- 
umn density along a given line-of-sight. Also gas in the disc, 
heated by feedback from super-nova explosions or photo- 
ionized, may contain O VII. Furthermore, satellite galaxies 
with ongoing star formation and related super-nova feed- 
back, as well as the photo-ionized surface layers of HI clouds 
in the halo (HVCs), may contribute to the O VII column 
density along some lines-of-sight. 

In our first set of N(OVII) estimates, in order to exclude 
such "contaminating" gas, we excluded all disc gas within a 
cylinder of radius r = Rd = 20 kpc, where Rd the galactic 
disc radius and height z — ±3 kpc, centred at the centre of 
the galaxy, and aligned with the disc. Moreover, it was found 
that potentially "contaminating" gas in satellites as well as 
in halo HI clouds was characterized by densities uh > nH,c, 
where nH,c = 10""^ cm~'^. On the other hand, the hot halo 
gas was found to everywhere have nn < riH.c- We hence 
excluded all gas in the halo of uh > nn.c from the analysis. 
The filling factor of this gas is very small, / <C 0.01, and the 
effect of excluding it on the results presented in the following 
is negligible. 

In order to assess the effect of removing the potentially 
contaminating gas we made a second set of N(OVII) esti- 
mates based on all gas, i.e., including gas in the disc, satellite 
galaxies and HI clouds. 

We should note here that, when choosing our lines of 
sight, as a first approach we did not try to achieve the same 
resolution in all areas of the celestial sphere. However, any 
over- or under- sampling of areas could result in differences 
in the calculated distributions. For this reason, we have com- 
pared between a constant 1 degree separation of the lines 
of sight and a simple weighting scheme that samples each 
galactic longitude with 360 • cos b lines, where b is the galac- 
tic latitude of the line of sight, and found little difference in 
the statistics. In the following we use the constant resolution 
for illustration purposes only and present all the analysis for 
the data calculated with the weighting scheme. 

Figure [2] shows the all sky, statistical N(OVII) distri- 
butions for the three galaxies and two types of estimates, 
employing the weighting scheme outlined above and figure 
[3] shows the cumulative histograms of the calculated column 
densities. 



Figure |4] shows contours of the logarithm of the O VII 
column density as a function of galactic longitude and lati- 
tude when all disc and "dense" (see above) halo gas particles 
are removed from the simulation data. Figure [5] shows the 
results with all the particles included. Both results have been 
calculated with a constant one degree distance between the 
lines of sight. 

For all galaxies, the halo only O VII column density 
is below the obs ervational upper limit of ~ 5 ■ 10^^ cm^^ 
(|Yao et al.ll200g ). There is a spread of approximately one 
order of magnitude in column density values that results 
from local fluctuations in density and/or metallicity. 

With the inclusion of all SPH particles the mean column 
density is not significantly affected, but the distributions all 
display a tail towards higher column densities (see Table 
[2 and figures [2] and [31) . The areas near the galactic plane 
are dominated by disc contribution, while halo objects give 
extended high column density areas at latitudes away from 
the galactic plane. 

As can be seen from figure O high O VII column den- 
sities are quite common when the disc is included. This is 
due to supernova explosions heating the gas and, of course, 
due to the high metallicity of the disc gas. It is worth not- 
ing that a feature near the disc plane of K26 can still be 
seen in figure |4l when the disc is removed. This is a typical 
example of the effects of supernova-driven outflow, driving 
high-metallicity hot gas into the lower halo. In the case of 
K33, the characteristic features in Figure[5]come from a large 
warp surrounding the galaxy, product of the tidal disruption 
of a satellite, previously accreted onto the galaxy. 

In order to assess how the inclusion of the cosmic UV 
background affects the results, we repeated the halo-only 
run using an oxygen fraction table calculated without a UV 
background. As illustrated in the left panel of Figure [21 the 
overall effect of the UV field is to shift the distribution to 
slightly lower column densities. This happens because the 
hardest photons of this radiation can ionise oxygen to levels 
higher than the sixth, so the latter remains relatively under- 
populated. However, this is not the case for K33, where the 
UV background has the opposite effect, that is to increase 
the column density values. This is easily understood, since 
the average temperature of the K33 halo is somewhat lower, 
due to its lower mass, so the UV field brings oxygen from 
lower ionisation stages to the sixth. 

The distance of 500kpc at which we decided to stop 
the integration was chosen so that we would only probe the 
halo of a single galaxy. Integrating to a somewhat smaller 
distance of 300kpc does not affect the results significantly, 
owing to the low gas densities at these distances from the 
galaxy. For distances larger than SOOkpc, one could consider 
the possibility that a line of sight passes through the haloes 
of neighbouring galaxies, which would in principle also give a 
contribution. To test the effect of a line of sight intersecting 
the halo of a neighbouring galaxy we calculated the N(OVII) 
values for 10000 random lines of sight passing through the 
halo of K15 (this galaxy was chosen because the virial tem- 
perature is comparable to that of M31 — see section 4). The 
lines of sight were chosen to originate and end at r=500 kpc 
distance from the central galaxy. Figure [6l shows the result 
of this calculation for impact parameters up to 300kpc. For 
larger impact parameters the contribution drops even fur- 
ther. 
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Table 2. Mean logarithm of column densities < log N(OV II) > 
for different runs. 
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It is apparent from Figure [6] that, unless a line of sight 
passes through a neighbouring halo with an impact parame- 
ter lower than ~100-200kpc, the contribution from this halo 
will be negligible. 



4 DISCUSSION AND CONCLUSIONS 

We have in this work determined the O VII column density 
of hot halo gas along various lines-of-sight in three simulated 
disc galaxies. In order to emulate observations performed at 
the solar position in the Milky Way, we have, in each galaxy, 
chosen a position in the midplane of disc, located 8 kpc from 
the centre of the galaxy. From this position we have drawn 
lines-of-sights emulating all sky coverage. 

We have tested that the results obtained do not in any 
significant way depend on where, along the R = 8 kpc mid- 
plane circle, the observer is assumed to be located. Moreover, 
we have for galaxy K33 tested that the results obtained do 
not in any significant way depend on the numerical resolu- 
tion of the galaxy formation simulation. 

For the two Milky Way sized galaxies, of Vc = 207 and 
245 km/sec, respectively, we find median halo gas O VII 
column densities of 1.5-2.5-10^* cm~'^, with a dispersion in 
log(A'^(0 VII)) of about 0.2. For the somewhat smaller disc 
galaxy K33, of Vc = 180 km/s, we find a median halo gas 
O VII column density of 3-10^^ cm~^, also with a disper- 
sion in \og{N{0 VII)) of about 0.2. The lower value for this 
galaxy is primarily due to the somewhat lower (virial) tem- 
perature of the hot gas (by about a factor of two), which, cf. 
Fig. [T] leads to a significant reduction of the number of oxy- 
gen ions in the sixth ionization stage for the relevant density 
range, nn ~ 10"^ - 10"^'^ cm"^. 

The main result of the paper is hence that for the 
two Milky Way sized galaxies, the median predicted halo 
O VII column density lies about a factor of 20 below the 
observational up per limit of 5 ■ 10^^ cm~^ reported by 
lYao et al.l (|2008h (and is comparable to the upper lim its of 
log{N{OVII)) < 14.2 - 14.8 set bv lYaoet alT(|201Gl ). but 
these limits are based on halo sight-lines, in general, pass- 
ing far away from the galactic centres a comparison is less 
straightforward — see also below). Moreover, not a single 
line-of-sight was found in any of t he galaxies with a halo 
O VII column density exceeding the lYao et al.l (|2008l ) limit. 

When the disc gas, as well as gas in satellites and HI 
clouds, is not excluded in the O VII column density esti- 
mates, the median A''(0 VII) increases by about a factor 
of two. Moreover, the A'^(0 VII) distributions display tails 
towards signi ficantly larger val ues with a few lines-of-sight 
exceeding the lYao et al.l (|2008l ) limit — this highlights the 
importance of excluding "contamination" from the disc and 



Figure 4. Logarithm of N(0 VII) in cm ^ for the three simu- 
lated galaxies, excluding the disc, satellites or HVCs. From top 
to bottom: K33, K26, K15 
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Figure 2. Histograms of ttie number of lines of sight at each column density bin. Left column: Blue dashed histograms result from 
including all gas particles, while red solid histograms result from excluding the disc and "dense" halo objects (see text for details). Right 
column: Red solid histograms are the same as shown in the right column (including a UV background and excluding the disc and "dense" 
halo objctcs), blue dashed histograms result from excluding the UV background (see text for details). Numbers are normalized to the 
total number of linos for each model. From top to bottom: K33, K26, K15 
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Figure 3. Cumulative histograms of the O VII column densities (fraction of lines of sight with column densities smaller than the 
corresponding column density bin). Left column: Blue dashed histograms result from including all gas particles, while red solid histograms 
result from excluding the disc and "dense" halo objects (sec text for details). Right column: Red solid histograms are the same as shown 
in the right column (including a UV background and excluding the disc and "dense" halo objetcs), blue dashed histograms result from 
excluding the UV background (sec text for details). From top to bottom: K33, K26, K15 
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Figure 5. Logarithm of N(0 VII) in cm -^ for the three simulated 
galaxies, including the disc and all objects in the halo. From top 
to bottom: K33, K26, K15 



Figure 6. N(0 VII) as a function of the impact parameter b for 
galaxy K15, which has a characteristic circular speed similar to 
that of M31. 
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other non-halo components, as discussed by the same au- 
thors. It also strongly suggests that the routinely reported 
ionic column densities of N{OVII) ~ lO^^cm"'^ at zero red- 
shift originate in gas located near or inside the galactic disc, 
while the extended gaseous corona has a much lower con- 
tribution. The largest values of 7V(0 VII) along a few lines- 
of-sight are mainly caused by either hot, super-nova driven 
bubbles in the disc gas or (in the case of K33) an outer, 
partly photo-ionized warp, originating from gas stripped off 
a previously accreted satellite. 

To test the sensitivity of our results on the assumed 
UVB and (simplified) radiative transfer scheme, we re- 
peated the calculations assuming no UVB. For the two Milky 
Way like galaxies, this leads to marginally larger median 
A'^(0 VII) values. The reason for this is that at the (virial) 
temperatures of the halo gas, O VII is the most abundant 
oxygen state, and with the UVB switched on, some of the 
O VII ions are photo-ionized to (mainly) the O VIII stage, 
lowering the population of O VII ions. For the lower mass 
galaxy K33, on the other hand, the effect of the UVB is to 
boost the O VII ion population. This is confirmed by calcu- 
lations (not presented in this paper) of the O VI and O VIII 
column densities for the same lines of sight. 

Finally, we have tried to quantify the effect a neighbour- 
ing galactic halo might have on the O VII column density 
for a random line of sight. The probability of that happen- 
ing and the relative effect this might have on the final re- 
sult depends, of course, on the sizes of the two galaxies and 
on the distance between them. In a simple calculation, we 
have randomly chosen the starting and ending coordinates 
of 10000 lines of sight outside a certain halo and calculated 
their O VII column densities. We find that the contribution 
is very small unless the impact parameter of the line of sight 
to the centre of the galaxy is lower than 100-200kpc. Includ- 
ing the disc gas and satellites in this case does not alter the 
results presented above unless the impact parameter is very 
small (smaller than 50kpc ) . Thi s is in agreement with what 
iBregman fc Llovd-Daviea ((20071) find for their line of sight 
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passing closest to M31, at a distance of 380kpc, giving no 
substantial contribution to the absorption. 

Considering this result, possible contamination to the 
lines of sight that probe the Milky Way halo could likely 
only come from M31 (more distant galaxies would have their 
O VII line redshifted, so it would be relatively easy to dis- 
tinguish their halo gas contribution from the one coming 
from local gas; this is beyond the scope of this paper, and 
will be the topic of a forthcoming paper — see also below). 
M31 is located at a distance of approximately 700kpc from 
our galaxy, which would cause it to contribute to the O VII 
column density with more than lO^'^cm"'^ for only about 
0.5% of the sky area (b < lOOkpc — see Fig. ^. Excluding 
b ^ 200kpc would reduce the M31 contribution to less than 
about 2 X lO^'^cm"^, and still only exclude about 2% of the 
sky area. 

As mentioned above, it will be the topic of a forthcom- 
ing paper to investigate the contribution from the haloes of 
other, more distant galaxies. This study will be based on 
the proper statistical descriptors at the large-scale galaxy 
distribution. 

In conclusion, the present observational upper limit 
is perfectly consistent with the results of state-of-the-art 
galaxy formation models, based on fully cosmological sim- 
ulations. Moreover, although the predicted halo A'^(0 VII) 
values lie below the observational upper limit, the difference 
is still only about an order of magnitude. With more sen- 
sitive future X-ray observatories, such as IXO, which will 
provide an order of magnitude increase in collecting area 
and a very high spectral resolution, it is very likely that one 
will be able detect O VII in the halo of the Milky Way, and 
possibly also in the haloes of other gala^cies. 
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